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CN . Abstract: Predictions from a new next-to- leading order (NLO) calculation for 

direct and resolved photoproduction of one and two jets are compared to simu- 
lated HERA data. We propose a method to match experimental and theoretical 
jet definitions and observe a reduced dependence on jet definitions and hadroniza- 
^ ■ tion corrections at larger transverse energies. From the irreducible uncertainty, we 

^ . estimate the maximum benefit that can be obtained from increased luminosity to 

\Q '. constrain the structure of the photon and the proton. 

^ 1 Introduction 

Oh! 

' Hard scattering of real photons off partons can be reliably predicted by perturbative QCD. 
The first NLO QCD calculation for complete (direct and resolved) photoproduction of one and 
^ two jets was completed recently It is based on the phase space slicing method and uses an 
^ , invariant mass cut to integrate soft and collinear singularities analytically. The program has 
■ ■ ' successfully been tested with an older existing program in single jet production 0. In order to 
be able to compare these theoretical cross sections to experimental data in photoproduction of 
jets from the ep collider HERA ^, the soft interaction links between initial and final state 
partons and hadrons have to be under control. Therefore it is only possible to extract new 
information on the parton densities in the proton and the photon if the jet definitions in the 
measurement and in the theoretical prediction match. In this paper we propose a method to 
achieve this for various jet algorithms. The experimental cross sections are simulated for 1994 
HERA conditions, where electrons of energy = 27.5 GeV collided with protons of energy 
Ep = 820 GeV, using HERWIG 5.8 [|. 

2 Jet Definitions 

According to the standardization of cone jet algorithms at the Snowmass meeting in 1990 0, 
calorimeter cells or partons i are included in a jet, if they have a distance of 

R^ = ^J{m - vjY + {<P^ - <PjY < R (1) 



from the jet center and a distance of 



from each other. Here, rj = — ln[tan(6'/2)] is the pseudorapidity related to the polar angle 
6, and is the azimuthal angle. If two partons have equal transverse energy they may be 
separated from each other by as much as 2R. As parton j does then not lie inside a cone of 
radius R around parton i and vice versa, one might with some justification also count the two 
partons separately. If one wants to study only the highest-i?T jet, this "double counting" must 
be excluded. The definition of an initiating cluster before a cone is introduced ( "seed- finding" ) 
is not fixed by the Snowmass convention, and different approaches are possible. The ZEUS 
collaboration at HERA uses two different cone algorithms: EUCELL takes the calorimeter 
cells in a window in 77 — space as seeds to find a cone with the highest Et- The cells in this 
cone are then removed, and the search is continued. On the other hand, PUCELL was adapted 
from CDF and starts with single calorimeter cells. It then iterates cones around each of them, 
until the set of enclosed cells is stable. In this case it may happen that two stable jets overlap. 
If the overlapping transverse energy amounts to a large fraction of the jets, they are merged, 
otherwise the overlapping energy is split. In addition, we simulate the same cross sections with 
the kT algorithm KTCLUS [0, §], where particles are combined if their distance 

d,,=mm{ET,,ET^fRl (3) 

is small. As the same recombination scheme is used, the results are quite similar to the PUCELL 
results. In the following we choose R = 1 throughout. Partonic jets with a large distance of 
two contributing partons are hard to find because of the missing seed in the jet center. This 
is especially true for the PUCELL algorithm, which does not perform a preclustering and does 
indeed find smaller cross sections and different hadronization corrections than the less affected 
EUCELL algorithm. We propose to model this theoretically by introducing an additional 
parameter i?sep to restrict the distance of two partons from each other [Q. This modifies eq. (j^) 
to 

Ex ~\~ Ex , 

^^Jp ^R^Rsep ■ (4) 

The meaningful range of i?sep is between 1 and 2. For two partons of similar or equal transverse 
energies Ex, Rsep is the limiting parameter, whereas it is the parton-jet distance R for two 
partons with large Ex imbalance. On a NLO three parton final state we find that Rsep = 1.5. ..2 
for EUCELL and R^ep ^ 1 for PUCELL and KTCLUS. 



Rij < min 



3 Results 

In figure |^, we show the dijet cross section da/dfj with Ex > 6 GeV and Arj = {rji — 772) G 
[—0.5, 0.5] for a) direct and b) resolved photoproduction as a function of the average pseudora- 
pidity of the two observed jets rJ = {rji +r]2)/2. The photon distribution in the electron is taken 
from the Weizsacker- Williams approximation with maximum virtuality of Q^ax — 4 GeV^ and 
longitudinal momentum fraction y G [0.2,0.8]. The direct and resolved regions are defined by 

x^^^ > 0.75 and x^^^ G [0.3, 0.75], where the sum in x'?^^ = ^'^^^ — runs over the two jets 




Figure 1: Demonstration of compatibility of the HERWIG Monte Carlo generator on matrix 
element level (histograms) and the leading-order (LO) perturbative Q CD prediction (curves) for 
a) direct and b) resolved dijet cross sections da/drj. 



with largest Et- In both regions, direct and resolved contributions are added because only 
their sum is physically meaningful. We demonstrate the compatibility of the two tools used 
in this study comparing the HERWIG Monte Carlo generator on LO matrix element level and 
the LO perturbative QCD prediction. We use the CTEQ3L proton and the GRV^(LO) photon 
structure functions. It was not possible to calculate as in 1-loop approximation in HERWIG, 
so we took the (inconsistent) choice of the 2-loop formula for with ^qqj) = 177 MeV for 

this part in the calculation as well. The scales in HERWIG could also not be changed from 
/i^ = = Mp = 2stu/ [s^ + + M^), but the effect with respect to using E^ as in the calcu- 
lation is very small. After these adjustments, HERWIG agrees with the LO QCD prediction. 

Figure ^ shows the same cross section for different jet algorithms. We compare our new 
NLO calculation with jet double counting, without jet double counting, and with Rsep = 1 
to simulated data from HERWIG with the EUCELL, PUCELL, and KTCLUS algorithms run 
on the final state particles. In the calculation, we now use a NLO set of input parameters, 
i.e. CTEQ3M proton, GRV^(HO) photon structure functions, and 2-loop as with ^qqj^ = 239 

MeV. For the Monte Carlo, we now take HERWIG including parton showers and hadronization. 
Due to the different parameters used in the NLO calculation and in HERWIG, we do not expect 
the overall normalization of HERWIG and NLO QCD to agree. However, the relative changes 
between no i?sep (equivalent to Rsep = 2) and i?sep = 1 on the theoretical side and EUCELL 
and PUCELL or KTCLUS on the experimental side show the expected similar behaviour, so 
that the i?sep parameter is well suited to match theoretical and experimental jet definitions. 
Jet double counting does not correspond to an experimental situation and is only shown to 
illustrate its effect on the theory. 

The E'r-dependence of the Snowmass jet definition uncertainties is shown in figure § for the 
a) direct and b) resolved dijet cross sections da/dEx, where we integrated over the complete 
rj range and over Ari G [—0.5,0.5]. The fractional difference of jet double counting from no 
double counting amounts to ~ 20% at 5 GeV and decreases continuously towards larger Et- 
Including the parameter Rg^p = 1 lowers the cross sections by as much as ~ 40% at 5 GeV, 
but its influence drops rapidly and gives a constant difference of about 20% (direct) and 10% 
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Figure 2: NLO (left) and HERWIG (right) predictions for direct (top) and resolved (bottom) 
dijet cross sections da /drj. We compare jet double counting (dotted), no Rgep (full), andRsep — 1 
(dashed) curves with EUCELL (full), PUCELL (dashed), and (dotted) histograms. 
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Figure 3: Ft -dependence of the Snowmass jet definition uncertainties for a) direct and b) 
resolved dijet cross sections da/dET- We demonstrate the effects of jet double counting (full 
curves) and setting i?sep = 1 (dashed curves) compared to no double counting with i?sep = 2. 



(resolved) almost over the whole i?T-range. Only at the boundary of phase space at very large 
the fractional difference increases again. Thus, one should match the jet definitions at 
small and large Et even more carefully. 
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Figure 4: Et- dependence of hadronization corrections for a) direct and b) resolved dijet cross 
sections da/dEx- 



The dependence of the hadronization corrections is shown in figure ^ for the a) direct 
and b) resolved dijet cross sections da/dET- We integrated again over the complete r/ range and 
over Arj e [—0.5, 0.5]. At low transverse energies of Et — 5 GeV, the hadronization corrections 
amount to ~ 20% - a comparable effect to the theoretical uncertainties discussed before. They 
decrease very nicely towards larger Et and vanish at 30 to 40 GeV to establish the correspon- 
dence between partonic and hadronic jets there. The error bars are due to limited statistics 
and could be drastically reduced with more computer time. 

Finally, we estimate the benefit that can be obtained from an increased HERA luminosity 
of 250 pb~^ to constrain the photon and the proton parton densities. Figure |^ gives the number 
of jets produced as a function of x^^^ (left), x^^^ (right), and transverse jet energy Et, where 
the bin sizes reflect the approximate expected experimental resolution. The cuts applied to 
the dijet cross section are y G [0.2,0.9], Et > 30 GeV, and rj < 2. If we require at least 
100 events, where statistical and systematic errors start to be of comparable size, jets with 
transverse energies up to 55 GeV can be measured, where the jet double counting uncertainties 
and hadronization corrections are very much reduced. We can still test the photon structure at 
large Et in the region of x^^^ = 0.4 — 1 and the proton structure in the region Xp^^ = 0.05 — 0.3. 



4 Conclusions 

Constraining the proton and photon structure functions in photoproduction of jets at HERA 
requires a good understanding of the jets in experiment and in theory. We used a new NLO 
calculation and simulated HERA data to match different experimental jet definitions (EUCELL, 
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Figure 5: Number of events produced at an increased HERA luminosity of 250 pb ^ as a function 
„o] 

S 



of x'^^^ (left), Xp^^ (right), and transverse jet energy Et- 



PUCELL, and KTCLUS) with theory predictions with different values for the i?sep parameter. 
At larger transverse energies, the uncertainties from different theoretical jet definitions and 
hadronization corrections are reduced. These regions can be studied if the HERA luminosity is 
increased to 250 pb~^, thus providing valuable information on the proton and photon structure 
functions over large x ranges. 
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